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1. Introduction 

In this paper we study the following full Magnet-Hydro-Dynamics system (MHD): 

' d t v + v ■ Vv — uAv + Vp = j x B 

d t E-V x B = -j 
(1.1) <d t B + VxE = 

divv = divB = 
K a(E + v x B) = j 

with the initial data 

v\ t =o = v , B\ t= o = B , E\ t =o = E . 

Here v, E, B: M.f x R 3 — > R 3 are vector fields defined on R 3 . The vector field 
v = (vi,V2,Vs) is the velocity of the fluid, v its viscosity and the scalar function p 
stands for the pressure. The vector fields E and B are the electric and magnetic 
fields of the fluid, respectively. And j is the electric current expressed by Omn's law. 
The force term j x B in the Navier-Stokes equations comes from Lorentz force under 
a quasi-neutrality assumption of the net charge carried by the fluid. Note that the 
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pressure p can be recovered from v and j x B via an explicit Calderon-Zygmund type 
operator. The second equation is the Ampere-Maxwell equation for an electric field 
E. The third equation is nothing but Farady's law. For a detailed introduction to the 
MHD, we refer to Davidson [2] and Biskamp pQ. 

Concerning the Cauchy problem associated to (11. II) . multiply the Navier-Stokes 
equations by v, the Ampere- Maxwell equations by (B, E) T and integrate (using the 
divergence free condition of the velocity) to get the following formal energy identity 

\j t [h\\l* + \\ B \\h + ll^liy + \\3\\h + "IIV* = o 

showing that both the viscosity and the electric resistivity dissipate energy. It also 
suggests that with initial data in (L 2 (IR d )) , one can expect to construct a global finite 
energy weak solution (a la Leray) . However, this intuitive expectation remains an 
interesting open problem for system (11.11) in both dimensions d = 2, 3 as an actual 
difficulty stands to derive some compactness especially for the magnetic field in virtue 
of the hyperbolicity of Maxwell's equations. This problem is similar to the global 
weak solvability of Euler's equations for ideal fluid which also remain open. In [9], 
P. L. -Lions suggested the notion of dissipative solutions for Euler's equation. These 
are functions which satisfy the energy inequality but not necessarily the system of 
equations itself. However, these solutions do coincide with any strong solution with 
the same initial data, if there is any. Dissipative solutions were also investigated for 
other MHD models as in [TJ]. For strong solutions, local wellposedness and small 
data global existence were investigated for other models of magneto-hydrodynamic as 
for example in [5] [6], but the nonlinearities there are weaker than the one in system 
1QJ) . In the two dimensional case with (v ,E ,B ) G L 2 (R 2 ) x (H S (R 2 )) 2 and s > 0, 
Masmoudi [12] proved the existence and uniqueness of global strong solutions to (11.11) . 
More recently, Keraani and the first author [7] showed the existence of global strong 
solutions in both dimension two and three with small initial data in spaces as close as 
possible to the energy space. However, the authors were not able to construct local 
solutions for arbitrary large initial data. 

Imposing more regularity on the initial electro-magnetic field, one can hope to solve 
(II. ip . In this paper we show the existence of local solutions with three dimensional 
large periodic initial data. Similar argument is valid in the case of the whole space and 
almost periodic case. We also show that if the initial data is not smooth enough, the 
solution is not smooth enough. In this consideration we use the periodic structure. 
Using the divergence free property of B, one can easily verify that 

V x (V x B) = -AB, 

and therefore, the magnetic field B satisfies an inhomogeneous damped wave equation. 
Let E and E be the divergence free and the gradient potential parts of the electric 
field E i.e. 



E := PE, and E := V(-A)~ 1 divE. 
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Note that E = E + E. In what follows, we use the following re-written equation 
coupled with damped wave equation: 

f d t v + P(v ■ Vv) - uAv = P(j x B) 
d tt B - AB + d t B = V x (v x B) 
(1.2) \ d a E - AE + d t E = -d t (v x B) 

d t E + E = -V(-A)^ 1 div(w x B) 
divf = divB = 



with the initial data 



*=o 



En. 



v\ t =o = v , B\t= = B , d t B\ t=0 = B 1 , E\ t=0 = PE , d t E\ t=0 := E\, E\ 

The formulation given by equation (II. 2p is good for constructing a mild solutions. Let 
C\ and £,2 be the propagators associated to the Fourier multiplier functions 



$i(*,ra) := e~* /2 cos ^\n\ 2 - 1/Atj , <5> 2 {t,n) := e" 



i/2_ 



sm 



VW - l/4t) 
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Mild solutions can be written as 



(1.3) u = u(sc,t) 



B = 5(x,t) 



e tA v 



e tA P [V • (v ® v) - E x B - (v x B) x B]ds 



e tA vn + M 1 ( y v,v) + M 2 (E,B)-M 3 (v,B,B) for x G T 3 , 
&{t)B + £ 2 (t)(B /2 + Bx) + [ C 2 (t - s) [V x (u x 5)] rfs 
C 1 (t)B + C 2 (t)(B /2 + B 1 ) + M 4 (v,B) for z G T 3 



E = E{x,t) = £ 1 {t)E + £ 2 {t)(E /2 + E 1 )- f £ 2 (t - s)d s P{v x B)(s)ds 

Jo 

= £ 1 {t)E + £ 2 {t){E /2 + E 1 ) + £ 2 {t){v xB )+ [ (d s £ 2 ){t - s)d s P(v x B)(s)ds 

Jo 

:= £ 1 (t)E + £2(t)(E /2 + E 1 ) + £ 2 (t)(v xB )-M 5 (v,B) for x G T 3 , 

and 

S = e-*E - / e- (< ~ s) V(A)- 1 div(t; x B) ds 
Jo 

:= E -M 6 (v,B) for x G T 3 , 

where I?i = <9 4 -B = —V x E . 

We handle periodic functions so that any smooth function can be decomposed as 

f{x) := f^'"'""- 

Throughout this paper, we assume that the mean value of the initial magnetic field is 
zero, namely Bq(0) = 0. The mean value of B is preserved along the original equation. 
Indeed, from 

d t B + V x E = 0, 

we have 

d t B(n,t) + m x E(n,t) = 
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which for n = gives dtB(0,t) = 0. If moreover 5(0,0) = 0, then B(0,t) = for any 
t > 0. 

In what follows, we consider the mild solution ( JTT3J) . To show local existence of (II. 3p . 
we need to find a function space in which one can show closed estimates for (II. 3p . In 
this point of view, L°° is one candidate of such function spaces since (jl.3p has a tri- 
linear term M 3 (w, B, B). More precisely, we can easily have the following estimate for 
M 3 : 

(1-4) ||M 3 |U < t\\{v xB)x BIU < tWvlUlBWl. 

However the L°° Lebesgue space is not suitable for the damped wave propagator C\ and 
£2- This difficulty is essentially the same as the unboundedness of singular integral 
operator in L°°. To overcome this, we introduce the weighted il £ 1 (in Fourier side) 
space with s weight given by 

X s := {u G S'(T 3 ) : \\u\\ X s = ||(1 + | ■ | 2 ) s/2 u||^ < 00}. 

We denote by := ||u||_yo. It is well known that X° is an algebra which is contin- 
uously embedded in BUC, the space of bounded uniformly continuous functions. For 
the periodic case, we can point out the following relationship between X s (T 3 ) and the 
Hdlder space C S (T 3 ). Let C S+0 (T 3 ) be C s+e (T 3 ) and X s+0 be X s+e for some e > 0. 

Proposition 1.1. We have u G C°°(T 3 ) if and only if u G r\ s >\X s . Moreover, 
C s+3 / 2+0 (T 3 ) is embedded in X s . 

The proof follows from the following well known inequalities, 

~ |M|ifo+3/2+0(T3) < ||u|| s s+3/2+0^ T 3^, 

where Bp q (T 3 ) is the Besov space. 

In this paper, we use the space X°. Since it is an algebra, one can easily estimate 
M 3 like (jl.4p . In addition, one can suitably estimate the damped wave operators L\ 
and £2 i n X Q . This kind of function spaces has already been used for fluid equations 
by several authors (see for example [21 IH El EH US])- However, in order to get a 
good representation of the electric field E as in (II ,2p . we need to decompose it into 
divergence free E and vector potential part E. Observe that we have 

ll-S'llxo < ||-E||x° + ll-E'llxo < ll-^llxo- 

The following local solvability and the propagator of C°° regularity is the first main 
result. 

Theorem 1.2. Let (v ,B ,E ) G (X° xl°x X°) with div v = div B = 0. Then 
there is a local-in-time unique mild solution of (11.11) 

(v,B,E) G C([0,T] :I°xI°x X°). 

Moreover, tf(v ,B ,E ) G C°° x C°° x C°° , then 

(v,B,E) G C°°([0,T] x R 3 ). 

Next we invest gate a "loss of smoothness" result. We can show that either the fluid 
flow or the magnetic field should develop a X 5+0 singularity within the lifetime of the 
solution even if initially neither the velocity field nor the magnetic field does have such 
singularities. This phenomena is caused by the damped-wave operator which does not 
have any smoothing effect (see [10] for example). We have the following second main 
result. 
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Theorem 1.3. For sufficiently small 5 > 0, any E 02 , -^0,3 £ C°° and (vq,B ) G 
C°° x C°°, there exist £ ,i G {X 1 \ X 1U ' 2 ) and a unique solution of (Oil on [0,T] 
such that one can choose time t G (0, T) with the property that either of the following 
alternatives can happen: 

v(t) £ X 5+<5 (T 3 ) or B(t) ^X 5+5 (T 3 ). 

Remark 1.4. For the usual 3D-Naveir-Stokes equations (which means B = E = 0), 
the solution is smooth enough even if the initial data vq is in X 1 \ X 1+5 ^ 2 . 

Now we give the rough idea to the proof of Theorem II .31 From the mild formulation 
of B with Bi — V x i?0; one can get inflation from the term LiB\. By choosing B 
smooth enough, CiB Q and £2-^0 remain smooth enough and if we assume that u{t) 
and B{t) remain smooth for time t, then ux B and £2(^)V(u x B) will be smooth and 
therefore can be absorbed by the inflation term. This leads to a contradiction. Thus, 
either term u(t) or Bit) must not be sufficiently smooth. 



2. Proof of the main results 
Proof of Theorem 1 1.2 . Let us set 

M t := max{ sup ||t>||, sup \\E\\, sup \\B\\} 

0<s<t 0<s<t 0<s<t 

Direct calculation shows the following estimates for Mj (j = 1, • • ■ 6): 



\\Mi(v,v)\\ 


< 


t 1/2 M? 


\\M 2 (E,B)\\ 


< 


tMl 


Ms(v,B,B)\\ 


< 


tMl 


\\M 4 (v,B)\\ 


< 


tMl 


11^5^,5)11 


< 


tMl 


\\M 6 (v,B)\\ 


< 


tM 2 . 



Indeed, on the one hand since 

sup s 1/2 |n|e" |n|2s < C, \\Pv\\ < C\\v\\, (v -V)v = V-(v®v ) and \\v®v\\ < \\v 

0<s<T 

we see that 

||Ml(«,«)|| < / SUp \\v{ S ')\\ 2 ds 

JO (J - s > ' 0<s'<T 

which clearly leads to the desired estimate of M\. Similarly, for M 2 , we see that 

\\M 2 (E,B)\\< [ sup \\E(s')\\\\B(s')\\ds, 

JO 0<s'<T 



6 



SLIM IBRAHIM AND TSUYOSHI YONEDA 



and this gives us the desired estimate for M 2 . The other terms are estimated similarly. 
On the other hand, we also have 



l|e* A n || 


< 


\\ u o\ 


£i(t)B \\ 


< 


\\Bo 


C 2 (t)B \\ 


< 


\\Bo 


C 2 (t)B 1 \\ 


< 

no 


\\Bi 


&{t)Eo\\ 


< 


\\Eo 


C 2 (t)E \\ 


< 


\\Bq 


£ 2 (t)£i|| 


< 


\\Ei 



By the above estimates, we can apply a usual fixed point argument. We sketch the 
proof in below and for full details refer for example to j3]. Define the successive 
approximation {%}j=i,2,- by 



ui(t) :-- 



B l 
Ex 
\E,J 



I 



e tA u 



C 1 (t)B Q + C 2 (t)(B /2 + B 1 ) 
A(t)£ + C 2 (t)(E /2 + E x ) + C 2 (t)(v x B ) 



V 



and 



u 



(t) :- 



( v j 
B 



-i\ 



j'+i 



'j'+i 



) 



(U:(c.r Vj) + M 2 (Ej, Bj) - U : .(c r B jt IS,)\ 
M 4 (vj, Bj) 
-M 5 ( Vj ,B 3 ) 
-MefaBj) 



Set 



and 



sup 

0<s<T 



sup ||-Ej(s) 

0<s<T 



sup 

0<s<T 



\Ejis) 



sup 

0<s<T 



\EM 



sup | 

0<s<T 



u 



j+1 



Using estimates (12. 1ft we have 

(2.2) Kj+i(T) <K X + CxT^K) + C 2 TK) + C 3 TKf. 

with positive constants Ci > 0, i — 1, 2, 3. The remaining part is now standard. From 
estimate (12. 2p . one can easily derive the uniform bound Kj(T) < 2K X < (||tt || + ||-B || + 
|| -Eb||) when T is sufficiently small with respect of the norm of the initial data. Similar 
calculation gives the pointwise estimate on Lj 

L j+1 (T) < CxT^KjLj + C 2 TK 3 L 3 + C 3 TK 2 Lj 

< 2C 1 T l/2 K 1 L j + 2C 2 TK 1 L j + 2C 3 TKfL j 

with positive constants Ci > 0, £ = 1,2,3. Now we take a sufficiently small T such 
that 2C 1 T 1 / 2 K l < 1/6, 2C 2 TK X < 1/6 and 2C 3 TK 2 < 1/6, we have 

L j+1 < (1/2)%. 

This gives us that there is a unique limit u such that Uj — > u G C([0, T] : X° x X° x X°) 

as j — > oo. It is easy to see that the limit it uniquely solves the equation. 

To show that (v, B, E) are smooth if (v , B , E ) are smooth enough, The argument 
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is, for example, similar to [TU Proposition 15.1]. From ( II. 3p and estimates (12. ip . we 
see that 

(2.3) ||Vv(t)|| < ||Vv || +^ (t-s)~ 1/2 ^||v||||Vv|| + \\E\\\\B\\ + ||u||||B|| 2 ^ ds 



\\VB{t)\\ < \\VB \\ + \\VB 1 \\ + J M|Vi;||||S|| + |M|||Vfl|M ds 
\\VE{t)\\ < \\VE \\ + \\VE 1 \\+ [ (\\Vv\\\\B\\ + \\v\\\\VB\\) ds 



\\VE(t)\\ < \\VE Q \\ +jf f||Vu||||£|| + IMII|Vfl||^ ds. 



By the above estimates (I2.3P and Gronwall's inequality, we can control || VJB|| + 1| Vu|| + 
|| V£|| + || V£|| uniformly in t G [0, T] since u, B,E,E G C([0, T] : X°). Repeating this 
argument for ||V fc -B|| + ||V fc t>|| + ||V fc J5|| + HV^II (k > 2), we can derive the space 
regularity of the solution. We can then get the time regularity from original equation 
(II. ip once we have the space regularity. □ 



Before proving Theorem ll.3| we need a pointwise estimate of M4. For that the 
following convolution estimate is the key. 

Lemma 2.1. For n G 1? , let p s {n) := C\j{C 2 + In^ 5 ), for some arbitrary positive 
constants C\ and C 2 . If s, s' > 0, then 

I p s (k)p s ,(n - k)\ < pmm{ s ,s>}(n) for all n G Z 3 . 
Proof. Since \n — k\ > \n\ — \k\ > \n\/2 for \k\ < \n\/2, we see that 

( Yl + E )\ps(k)pAn-k)\ < J2\p s (k)p s ,(\n\/2)\ + J2\Ps(\ n \WPs'( n 

fcez 3 ,|fc|<|n|/2 fcez 3 ,|fc|>|n|/2 fcez 3 fcez 3 

< m&x{p s (n),p sl (n)}. 

□ 

Remark 2.2. We see that if u G" X s ~°, then there is {n.,}^ — > 00) such that 
|^(^j)| > Ps{ n j)- Moreover, if u G X s , then |w(n)| < p s - 3 (n) f° r a h fi G Z 3 . 

Proof of Theorem I Recall that 

(n 3 E 0t2 (n) - n 2 E 0i3 (n) 
ni E 0t3 (n) - n 3 E 0A (n) 
n 2 E ,i(n) - nxEo^in) 

We focus on the second component of the vector in (12 .4p . We choose E^i in order to 
satisfy |-Eo,i(h)| > Pi+s/ 2 {n). In particular, this gives us that \n 3 E 0>1 (n)\ > pi+8/ 2 (n) 
for all n 3 7^ 0. 

Assume that (v,B) G C([0,T] : X 5+(5 (T 3 )). We show that a loss of regularity is 
brought from non-smoothness of Eq not Bq. By Remark 12.21 and Proposition 11.11 we 
have 

\v(t, n)\, \B(t, n)\ < p 2 +<5-o( n ) with eventually C\ = C\{T) and C 2 = C 2 {T). 



By Lemma [2.11 we see that 

\F [v x B] (t, n) I < p 2+ s-o{n) for t < T. 
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thus, 

\M A {v,B)\ < 7> 2+<5 _ (n). 
On the other hand, for any n (with n 3 7^ 3), there is t G (0,T) (in this case we set 
t = ( 2 ^ +1 )/ 2 7r < y f or S ome fceZ) such that 

I $2(^,^)^3^0,1 (n)| > p2+s/2{n). 

Since So, -^0,2 and £0,3 are smooth enough, we see that 

|$i(i,70(£oW/2)|, (^M^oWl, |$ 2 (t, 71)71^0,3(71) I < p a (n) 

for any s > (in this case we set s = 2 + 5 — 0) and t G [0, T]. Thus for any n (with 
n 3 ^ 0), there is t G (0, T), we have the following lower bound for the second vector of 
B: 

\B 2 (t,n)\ > \§ 2 (t,n)n 3 E 0A (n) \ - p 2+s _ (n) -Tp 2+S _ (n) > p2+s/2(n). 

This is a contradiction against B{t) G C([0,T] : X 5+<5 (T 3 )), since |JB(t, n)| must be 
controlled by p 2 +<5-o( ri ) f° r an Y t G [0, T]. □ 
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